Micro-slip field at a rough contact driven towards macroscopic sliding 



J. Scheibcrt, G. Dcbregeas, A. PrcvoslQ 
Laboratoire de Physique Statistique de I'ENS, UMR 8550, 
CNRS-ENS-Universites Paris 6 and 7, 24 rue Lhomond F-75231 Paris, France 
(Dated: September 18, 2008) 

The incipient sliding ol a multicontact interface between a smooth glass sphere and a rough 
elastomer block is studied using an image correlation velocimetry technique. The displacement field 
of the elastomer's surface is measured both inside and outside the apparent contact region with a 
submicrometer displacement resolution. For a given shear load, a coexistence between an inner stick 
region surrounded by an outer slip annulus is observed. The stick region extension decreases with 
increasing load and eventually vanishes at the sliding threshold. These measurements allow for the 
first quantitative test of Cattaneo and Mindlin (CM) classical model of the incipient sliding of a 
smooth interface. Small deviations are observed and interpreted as a result of the finite compliance 
of the rough interface, a behavior which contrasts with Amontons' law of friction assumed to be 
valid locally in CM's model. We illustrate how these measurements actually provide a method for 
probing the rheology of the rough interface, which we find to be of the elasto-plastic type. 
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Amontons' law of friction states that the shear force 
required to trigger relative sliding between two solids in 
contact is proportional to the applied load normal to the 
interface. However, it is known that even for minute 
shear forces, well before macroscopic sliding settles in, 
micro-slip occurs within the contact. Slip bearing regions 
are thus expected to coexist with adhesive ones and to 
progressively invade the contact as the shear force ap- 
proaches the sliding threshold. Providing a proper local 
description of this incipient sliding regime remains chal- 
lenging despite its importance to the fields of tribology 
1 1 , earthquakes nucleation Q or mechanical engineering 
3, In the case of spherical contacts between elas- 
tic bodies (Hertz contact), an analytical solution to this 
problem is given by Cattaneo and Mindlin's (CM) pro- 
cedure which has now become a classical model for shear 
loaded frictional contacts below the sliding threshold. It 
predicts that the stick region is a central disk whose ra- 
dius continuously decreases from the apparent contact 
zone radius down to zero upon increasing shear loading 
0, @|- Experimental confirmation of CM's predictions 
has been obtained only indirectly by measuring macro- 
scopic force-displacement curves or by observing the fret- 
ting zone resulting from the interfacial micro-slip Q . 

In CM's approach, the interface is assumed to be 
smooth, incompressible and to obey locally Amontons' 
rigid-plastic like friction law. This is at odds with most 
realistic situations where surfaces are rough down to the 
micrometer scale, resulting in a multicontact interface 
with a much richer mechanical response [1, [^, [l^, [ll|, [l3] ■ 
One may wonder how these characteristics modify CM's 
predictions. In this Letter, we address this question by 
probing the local displacement field at a contact between 
a rigid sphere and a rough elastomer block for which dis- 
placements can be measured with a submicrometer reso- 



lution using an image correlation velocimetry technique. 
This method directly makes use of the micro-junctions 
and micro-asperities present at the interface as markers. 

The sphere is an optically smooth glass lens (radius 
of curvature R = 128.8 mm) and the elastomer block 
(50 X 50 mm, thickness h ^ 15 mm) is made of a 
crosslinkcd PolyDimethylSiloxane (PDMS Sylgard 184, 
Dow Corning, Young's modulus E = 2.2± 0.1 MPa, Pois- 
son ratio i/=0.5). Its surface has been rendered rough by 
moulding the cross-linker /PDMS melt against a Plexiglas 
surface mechanically abraded with an aqueous solution of 
SiC powder. Its topography was characterized with op- 
tical profilometry (M3D, Fogale Nanotech). Its power 
spectral density is a power law down to the micrometer 
scale and up to about 40 /im. Its characteristic thickness 
(T; taken as the standard deviation of the height distri- 
bution, was found to be 1.28 ± 0.05 fim. Both glass and 
PDMS surfaces were cleaned with ethanol and dried with 
filtered air prior to any experiment. 

A typical experiment consists in pulling the elastomer 
block in contact with the sphere under a prescribed nor- 
mal load P at constant velocity V with a precision motor- 
ized actuator (LTA-HL, Newport), while recording simul- 
taneously P, the shear load Q and images of the interface 
(Fig. [T|). In all experiments, V ~ 4 /im/s, a velocity small 
enough for visco-elastic interfacial dissipation to be neg- 
ligible [l3|. Both P and Q are measured at a sampling 
rate of 1 kHz by probing with position sensors the deflec- 
tions of two cantilevers (Fig. [1] lower inset). The normal 
loading of the contact produces a significant shear force 
due to the coupling between normal and lateral motion 
of each cantilever. To circumvent this problem the con- 
tact is manually renewed until Q is less than 1% of P. 
During this separation procedure, no measurable pull-off 
force is observed, indicating that adhesion forces are neg- 
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FIG. 1: (Color online) Sketch of the experimental setup (lower 
inset), time evolution of the shear force Q during the loading 
experiment with P = 0.33 N and V = i fxm/s, and image 
of the contact (upper inset) with its apparent contact border 
(white dashed line) and coordinates system. On the sketch, 

(1) is the shear cantilever (stiffness 9579 ± 25 N.m~^) and 

(2) the normal one (stiffness 689 ± 5 N.m~^), (3) and (4) are 
two capacitive position sensors (MCCIO and MCC20, Fogale 
Nanotech), (5) is a glass plate to which the PDMS block (6) 
is attached, (7) is the glass lens and (8) the resulting circular 
contact. On the Q curve, circled points and black disks indi- 
cate Q values at which the displacements are displayed in the 
next figures. 



ligible [ij, ll5[. Q versus time curves, like the one shown 
in Fig. [T] for P = 0.33 N are reproducible for any P in 
the experimentally accessible range [0, 1 N], with the fol- 
lowing typical behavior. While P varies by less than 1%, 
Q increases until it reaches a maximum value Qs beyond 
which it slightly decreases before flattening at a constant 
value Qd = 0.96 Qs signaling a steady sliding regime with 
a dynamical coefficient of friction of about 0.6. 

Imaging of the contact is done by illuminating the op- 
tically transparent PDMS block from below with a white 
LED and a ground glass diffuser, and using a stereomicro- 
scope (Olympus SZll). Images (Fig. [TJ upper inset) are 
recorded by a CCD camera (Hamamatsu C8484-05G, 12 
bits), whose 1344 x 1024 pixels sensor provides at the cho- 
sen magnification a field of view of 9.8 x 7.5 mm. Images' 
contrast results from the diffusive nature of the rough 
interface. In the contact region, additional bright spots 
correspond to the micro-junctions which favorably trans- 
mit light. Intcrfacial displacement fields are extracted 
from snapshots acquired at a 4 Hz frame rate using an 
image correlation velocimetry technique [1^, 17 1 . 

It consists in finding, for a given subimagc at posi- 
tion (x, y) in a reference image the displacement (u^, Uy) 
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FIG. 2: (Color online) Snapshots of the displacement field 
Ux{x,y) at the positions labeled 0, 1, 2 and 3 on the curve 
of Fig. [T] The reference, undeformed image is taken at posi- 
tion 0. u{x,y) is represented by joining together with a solid 
line the extremities of the displacement vectors Ux of equal 
abscissa grid points. Since lies in the range [0, 100 /im] an 
amplification factor of 30 has been used to ease vizualisation. 
The dashed circle is the border of the apparent contact region 
whose area is unaffected by the loading. The grey circular disk 
is a region of null displacement within the contact. 



which provides the maximum correlation with a subse- 
quent deformed image. Since Ux is at least an order of 
magnitude larger than Uy, the correlation function is cal- 
culated by translating a square window of side length A 
along the shear direction x only. With A = 20 pixels, 
the Ux field is obtained with a spatial resolution of about 
150 ^m. The resolution on the displacements, found by 
correlating two different images of the interface prior to 
contact, and taken as the standard deviation of the dis- 
placements distribution, is about 275 nm, i.e. 1/30*'' 
pixel. All displacements are defined respective to an im- 
age taken at Q = N. 

The displacement field at various applied forces along 
the loading curve of Fig. [l]was determined by correlating 
successive images with the Q = Q reference image and is 
shown in Fig. [2] at four instants. For a given Q a central 
circular stick region coexists with a surrounding annulus 
in which slip has already occured. Its diameter decreases 
with Q and vanishes at Qs- When Q = Qd, Ux is the sum 
of two terms. The first one corresponds to the driving 
motion and is a uniform displacement of amplitude Vt. 
The second one is a residual displacement corresponding 
to the deformation of the elastomer block (Fig. [31 insets). 

Such a scenario has been theoretically predicted by 
Cattaneo and Mindhn (CM) 0, Q who have considered 
the combined normal and shear loading of a contact be- 
tween spherical elastic bodies. Their calculations assume 
that (1) both surfaces are smooth, (2) the pressure distri- 
bution p within the contact is unchanged upon shearing 
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and given by Hertz contact theory, and (3) Amontons' 
law of friction is valid locally at any (r, 9) position (Fig. 
m upper inset), i.e. slip occurs wherever the shear stress 
q reaches fip, /i being the friction coefficient. For a con- 
tact of radius a between a rigid sphere of radius R and an 
elastic plane of reduced modulus E* = under nor- 

mal load P, p{r, 6) is given by p{r) ~ po\/l — j with 
P% = CM predicts the coexistence of an inner 
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adhesive circular region of radius c = a (1 — Q/(^P)) 
surrounded by an outer slip annulus. Using a superpo- 
sition principle, CM's calculations provide complete an- 
alytic expressions for Ux within the contact, in both the 
stick and slip regions 0, 01 • Further derivations by John- 
son 0] also give Ux outside the contact, thus providing the 
entire field at the interface. Comparison between CM's 
and measured Ux were done in two ways, by averaging 
Ux{r,9) over 9 (Fig. [31[a)), and by evaluating Ux{a,9) 
with 9 running from to 27r (Fig. [3)Jb)). 

In Fig. [3ja), displacement curves are shown for all 
values of Q corresponding to the dotted positions along 
the Q curve in Fig. [1] The inset shows Ux{r) — Ux{0) in 
the steady sliding regime at the last four dotted points. 
Evaluation of Ux{a,9) (Fig. [3Kb)) is done similarly at 
the same dotted points in the transient loading regime 
and in steady sliding (inset). As shown, a good over- 
all agreement between CM's predictions (red solid lines) 
and the measured Ux is found. Close look at their radial 
dependence even shows that the radius of the stick re- 
gion is very close to CM's prediction for c. In addition, 
CM's expressions for Ux{a,6) predict a 005(26*) depen- 
dence, which is well reproduced over the whole 9 range 
and for all Q. 

Systematic deviations to CM's predictions are how- 
ever clearly seen on the radial profiles displayed in Fig. 
[Sja). Outside the stick region (r > c), the experimental 
data points systematically lay below the predicted pro- 
files. Furthermore, while CM's model predicts a kink in 
the radial displacement profiles at r = a, the measured 
curves remain smoother. This may be attributed to de- 
viations to Hertz's pressure field around r = a, induced 
by the non-linear normal compressibility of the rough 
layer More surprisingly, in the stick region (r < c) 
a non-zero displacement in the submicromctric range is 
measured whose value increases with Q. This effect can 
be best caught by evaluating Ux{r = 0) averaged over a 
disk of radius A centered on the contact as a function of 
time. Results plotted in Fig. HJa) and magnified on Fig. 
HKb) show that in the loading phase Ux{r = 0) increases 
continuously from to about 1 fim. 

This can be understood when considering that the im- 
age correlation velocimetry technique actually probes the 
mean displacement over the mean thickness a of the 
rough layer. It indeed averages out intensity fluctuations 
due to both the micro-contacts and the non-contacting 
micro-asperities between micro-contacts. We can then 
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FIG. 3: (Color online) Comparison between the measured 
displacements (solid lines with filled disks) and CM's predic- 
tions (red solid lines), (a) Ux{r)/a versus r/a. All curves are 
shifted arbitrarily along the y axis for vizualization and CM's 
curves have been smoothed out with a running window of size 
A. From bottom to top, each set of curves corresponds to an 
increasing Q in the loading phase (black dots in Fig. [!}. In- 
set shows 4 almost fully overlapped Ux{r)/a curves in steady 
sliding measured at the last 4 dotted positions in Fig. [T] (b) 
Ux{a,8)/a versus 8 evaluated at the same Q in the transient 
regime, and in steady sliding (inset). 



propose the following scenario. As long as Q <Qs, sum- 
mits of the asperities in contact remain stuck [22l | while 
their underlying bulk base is displaced due to the ap- 
plied shear. The rough layer is therefore deformed with 
a mean strain which can be estimated, at the center of 
the contact, as eq = Ux{r = 0)/a (Fig. HKc)). The image 
correlation velocimetry technique developped here thus 
appears as a tool to probe locally the rheology of the 
rough layer, provided that the corresponding shear stress 
go is known. Since our experiment does not allow any 
direct shear stress measurement, we have used CM's pre- 
diction for go = qir = 0) = /ipo(l — c/a) to obtain go 
versus eo (Fig. HKd)). Clearly, the rough layer mechan- 
ical response deviates significantly from the Amontons' 
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the classical model of Cattaneo and Mindlin thus pro- 
viding its first direct experimental test. Significant de- 
viations are however observed and shown to result from 
the too stringeant smoothness assumption, which is of- 
ten made in continuum mechanics calculations. The in- 
terfacial rough layer indeed exhibits a finite shear com- 
pliance. The latter was directly probed, allowing to ev- 
idence an elasto-plastic like friction law different from 
Amontons' law which is commonly believed to apply for 
multicontacts. These results provide direct experimen- 
tal support to the interpretation of a recently observed 
non- Amontons behavior in a similar setup 21 1 . 
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FIG. 4: (a) u^(r = 0) versus time. The dashed line is a 
linear fit, characterizing steady sliding at V = 4 /im/s. (b) 
Magnified centered region of the previous plot, (c) Sketch 
detaining the mechanism resulting in a non-zero measured 
displacement in the stick region, (d) Calculated shear stress 
go at r = using CM's model versus the strain eo defined 
as Ux (r = 0)/(7. Shown in dashed lines are the go ~ 20 kPa 
horizontal line and a linear fit of the data points for eo < 1- 



rigid-plastic like behavior assumed by CM. The friction 
law is rather elasto-plastic like with an effective shear 
modulus Gi given by G; = dqQ/2deQ of about 10 kPa. 
This value is two orders of magnitude smaller than the 
bulk shear modulus G = 2(1+^) « 0.73 MPa. It is fully 
consistent with previous measurements on multicontacts 
averaged over the whole apparent contact area ■ It is 
expected to depend on the local pressure so that system- 
atic measurements at successive different loads would be 
required to determine the complete mechanical response 
of the rough layer. 

In this Letter we have implemented an image correla- 
tion velocimetry technique to measure with a submicrom- 
eter resolution the displacement field at a multicontact 
interface between a soft elastomer and a rigid body. This 
technique provides the first direct and non-invasive mea- 
surements of the interfacial slip field in both incipient 
sliding and steady sliding regimes at a sphere-on-plane 
contact. It can be directly applied to dynamical regimes, 
therefore complementing real contact area measurements 
on rough contacts such as those recently implemented 
in Plexigias to investigate the interfacial dynamics of a 
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multicontact interface at the onset of sliding 
principle, it could also be extended to smooth contacts 
by patterning the elastomer with markers located below 
its surface [20| 

The measured displacement fields are well captured by 
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